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ON . 

Factorization scheme analysis of ^(x, Q 2 ) in the next-to-leading order QCD is revis- 
ited. It is emphasized that the presence of the inhomogeneous term in the evolution 
Q_e equations for quark distribution functions of the photon implies subtle but important 

difference in the way factorization mechanism works in photon-hadron and photon- 
■ photon collisions as compared to the hadronic ones. It is argued that none of the 

t-h . existing NLO analyses of F^x^Q 2 ) takes this difference properly into account. The 

source of the ensuing incompleteness is traced back to the misinterpretation of the be- 
haviour of q' y (x,M) as a function of a s (M). Parton model interpretation of the so 

called "constant terms" in the LO photonic coefficient function C^°\x) is given and 
— |- . smooth transition between the properties of virtual and real photon analyzed. Finally 

t-h I phenomenological consequences of this analysis are discussed. 

00 ■ 

1 Introduction 

^ : 

Observed from a large distance the photon behaves as a neutral structureless object governed by the 
laws of Quantum Electrodynamics. However, when probed at short distances it exhibits also some 
properties characteristic of hadrons. This "photon structure" is quantified, similarly as in the case 
of hadrons, in terms of parton distribution functions (PDF), satisfying certain evolution equations. 
Because of the direct coupling photons to quark-antiquark pairs these evolution equations are, 
contrary the case of hadrons, inhomogeneous. This inhomogeneity has important implications for 
the way factorization of mass singularities operates in collisions involving photons, implications that 
have not been properly taken into account in existing NLO analyses of i^(x,Q 2 ). The primary 
aim of this paper is to remove this shortcoming. 

Secondly, we shall address several issues concerning the structure of the virtual photon: tran- 
sition between the properties of real and virtual photon, properties and role of the longitudinal 
virtual photon, and parton model interpretation of the so called "constant terms" in LO photonic 
coefficient function C^\x). 

The paper is organized as follows. In the next Section basic facts and notation concerning PDF 
of the photon are reviewed and the properties of the pointlike part of quark distribution function 
of the photon critically reanalyzed. In Section 3, which contains the main result of this paper, we 
discuss in detail the factorization scale and scheme dependence of F^(x, Q 2 ) at the NLO and point 
out the ingredients that must be included to make this analysis complete. The properties of the 
virtual photon and the transition of PDF of the virtual photon to those of the real one are analyzed 
in Section 4. Phenomenological implications of the present analysis are discussed in Section 5. 
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2 Structure of the real photon 



Despite the recent progress in investigation of the structure of the photon 2 our knowledge of the 
properties of the photon still lags behind that of the nucleon. We shall be primarily interested 
in strong interaction effects, but as the basic ideas and formalism of the partonic structure of the 
photon have a close analogy in QED, the latter will serve as a guide in some of the following 
considerations. 



2.1 Notation and basic facts 

In QCD the coupling of quarks and gluons is characterized by the renormalized colour coupling 
("couplant" for short) a s (fi), depending on the renormalization scale \x and satisfying the equation 

da s (n) p 2 Pi 3 

j^jj = 0(a t (p)) = -^<* a (ji) - Y^^) + ■■■, (1) 

where, in QCD with rif massless quark flavours, the first two coefficients, (3q = 11 — 2nj/3 and 
Pi = 102 — 38nj/3, are unique, while all the higher order ones are ambiguous. As we shall stay in 
this paper within the NLO, only the first two, unique, terms in ([!]) will be taken into account in the 
following. Nevertheless, even for a given r.h.s. of (|]) its solution a s (/u,) is not a unique function of 
fi, because there is an infinite number of solutions of (|l]), differing by the initial condition. This so 
called renormalization scheme (RS) ambiguity f\ can be parameterized in a number of ways. One of 
them makes use of the fact that in the process of renormalization another dimensional parameter, 
denoted usually A, inevitably appears in the theory. This parameter depends on the RS and at the 
NLO actually fully specifies it: RS={Ars}. For instance, a s (/i) in the familiar MS and MS RS are 
solutions of the same equation (|l|), but are associated with different Ars Q- In this paper we shall 
work in the standard MS RS of the couplant. 

In QCD "dressed" PDF f\ result from the resummation of multiple parton emissions off the cor- 
responding "bare" parton distributions. As a result of this resummation PDF acquire dependence 
on the factorization scale M. In parton model this scale defines the upper limit on some measure 
t of the off-shellness of partons included in the definition of D(x, M) 

Di(x,M)= dt(k(x,t), i = q,q,G, (2) 



where the unintegrated PDF di(x, t) describe distribution functions of partons with the momentum 
fraction x and fixed off-shellness t. Parton virtuality r =| p 2 —m 2 | or transverse mass m\ = p^+m 2 , 
are two standard choices of such a measure. Because at small t, di(x,t) = 0(l/t k ),k = 1,2, the 
dominant part of the integral (||) comes from the region of small off-shellness t. Varying the upper 
bound M 2 in (|2|) has therefore only a small effect on the integral @, leading to weak (at most 
logarithmic) scaling violations. 

The factorization scale dependence of PDF of the photon f] is determined by a system of coupled 
inhomogeneous evolution equations 

dS(x,M) 



d In Ml 



k q + P qq ®Y,+ P qG (8) G, (3) 



2 For recent theoretical and experimental reviews see |lj and ||, respectively. 
3 In higher orders this ambiguity includes also the arbitrariness of the coefficients /3i, i > 2 in dl]). 
4 The variation of both the renormalization scale fj, and the renormalization scheme RS={Ars} is actually redun- 
dant. It suffices to fix one of them and vary the other, but 1 will stick to the common habit of considering both of 
them as free parameters. 

5 In the following the adjective "dressed" will be dropped. 

6 If not stated otherwise all distribution functions in the following concern the photon. 
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dG(x, M) 
dlnM 2 
dg^ s (x,M) 



k G + Pa q 8> S + P GG ® G, 

(JNsfcg + P N S ® 9NS ) 



dlnM 2 

where cr^g = (e 2 /(e 2 ) — l)/rif and the convolution (8) is defined in a standard way as 

- 1 dy 



P <8> q{x) 



x y 



-P(x/y)q{y). 



(4) 
(5) 

(6) 



The singlet and nonsinglet quark distribution functions S and qf are given as 

n f n f 



i=l 



8=1 



^ s (x,M) = qt(x,M 2 )-E(x,M)/n f . 



(7) 
(8) 



To order a the splitting functions Pij(x, M) and M) are given as power expansions in a s (M): 



a 



k q (x,M) = — 



k G (x,M) 
Pij{x,M) 



a 
2^ 



«f(«) + ^*fw + (^)Vw+- 

^g) W+ (-M)%g) W+ ..: 



.(0) 



2tt 



0) 
(10) 

(11) 



where the leading order splitting functions k q v> (x) = 3e 2 (x 2 + (1 — x) 2 ) and P^\x) are unique, 

while all higher order ones kg , Kq , Pu, , j > 1 depend on the choice of the factorization scheme 
(FS) []. The photon structure function F$(x, Q 2 ), measured in deep inelastic scattering experiments 
on the (slightly off-shell) photons || || || , is given as the convolution f\ 



F 2 7 (x, Q 2 ) = J2 ^e 2 \q{M) ® C q (Q/M) + G[M) ® C G {QJM) + C*(Q/M) 



(12) 



of photonic PDF and coefficient functions C q (x), C G (x), C!*(x) admitting perturbative expansions 



C q {x,Q/M) 
C G (x,Q/M) 
C«{x,Q/M) 



5(1 -x) 



+ 



2?r 



d 1 )(x,g/M) + -.., 



2tt 



C^(x,Q/M) + - 



)^\x,Q/M) + e 2 q ^^C?( X ,Q/M) 



+ 



(13) 
(14) 
(15) 



The renormalization scale /U, used as argument of the expansion parameter a s (/x), is in principle 
independent of the factorization scale M. Note that despite the presence of /i as argument of 
a s (fi) in (1~3|-15)) the coefficient functions C q , C G and are actually independent of \i because the 

/i-dependence of a s {ji) is cancelled by explicit dependence of C q , Cq , Cy\i > 2 on \i Q. On the 

7 We can turn this statement around and consider any factorization scheme FS to be specified by the corresponding 
set of functions k ( q j) , , pff , j > 1. 

8 The factor 2 accounts for the inclusion of antiquarks in the sum. 
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other hand, PDF and the coefficient functions C q , Cq and do depend on both the factorization 

scale M and factorization scheme FS={kq , Kq , , i > 1}, but in such a correlated manner that 
physical quantities, like i 7 ^, are independent of both M and the FS, provided expansions (|9 11 ) and 



(13|-15|) are taken to all orders in a s (M) and a s (/x). In practical calculations based on truncated 
forms of (^-11) and (p^-15) this invariance is, however, lost and the choice of both M and FS makes 
numerical difference even for physical quantities. At the NLO RS={A RS } and FSHX^ k ? \ P i?}- 

The expressions for C { q\c ( (P given in Q are usually claimed to correspond to "MS factorization 
scheme". As argued in ||, this denomination is, however, incomplete. The adjective "MS" concerns 
exclusively the choice of the RS of the couplant a s and has nothing to do with the choice of the 
splitting functions p[j ■ The choices of the renormalization scheme of the couplant a s (M) and of 
the factorization scheme of PDF are two completely independent decisions, concerning two different 
and in general unrelated redefinition procedures. Both are necessary in order to specify uniquely 
the results of fixed order perturbative calculations, but we may combine any choice of the RS of 
the couplant with any choice of the FS of PDF. Note that the coefficient functions Cq \ Cq , C*p 
depend on both of them, while the splitting functions depend only on the FS of PDF. The results 
given in § correspond to MS RS of the couplant but to the "minimal subtraction" FS of PDF |]. It 
is therefore more appropriate to call this full specification of the renormalization and factorization 
schemes as "MS + MS scheme" . Although the phenomenological relevance of treating \x and M 
as independent parameters has been demonstrated (T^j, we shall follow the usual practice and set 
/i = M. 

2.2 Properties of the pointlike solutions 



The general solution of the evolution equations (|3]-|5|) for the quark distribution functions q(x, M) 
can be written as the sum of a particular solution of the full inhomogeneous equation, called pointlike 
part, and the general solution of the corresponding homogeneous equation, called hadronic part: 

q(x,M)=q PL {x,M)+q UAD (x,M). (16) 

Having separated q(x, M) into its hadronic and pointlike parts we can insert it into the r.h.s. of the 
evolution equation (||) for the gluon distribution function and write its solutions as a sum of the 
hadronic and pointlike parts. At the LO, where kc = 0, both of these components satisfy standard 
homogeneous evolution equation 

^P = P G ,®E i + P GG ®G i , k = HAD, PL (17) 

However, as there is an infinite number of pointlike solutions q PL (x, M), which differ by terms 
satisfying the homogeneous evolution equation, the above separation of the quark and gluons dis- 
tribution functions into their pointlike and hadronic parts is not unique and consequently these 
concepts have separately no physical meaning. To see the most important feature of the pointlike 
part of quark distribution functions that will be crucial for the following analysis, we now consider 
in detail the case of nonsinglet quark distribution function q^ s (x, M) at the LO (which is sufficient 
for our purposes) and, moreover, drop the superscript i. A subset of all pointlike solutions of (||) 
may be characterized by the value of the initial scale Mq at which they vanish. At the LO and in 
terms of moments such solutions are given explicitly as 



q^(n,M ,M) 



a s {M) 

? See Section 2.6 of |ic|, in particular eq. (2.31), for discussion of this point 



as (M) y 2P ™^ 



a NS (n), (18) 
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Figure 1: Diagrams defining pointlike parts of nonsinglet quark and gluon distribution functions of 
the photon in LL approximation. 



where 



a in 



a NS (n) = -— . (19) 

and result from resummation of infinite series of diagrams in the upper part of Fig. |]: 

which, as illustrated in Fig. 2, softens substantially the x— dependence of a^s(x) with respect 
to the first term in (pQ), proportional to /cnsW- This construction is similar to that of hadrons 
with the first term in (|20|) , (a/2iT)k^(x) Iyi(M 2 /Mq), playing the role of "bare" quark distribution 
function. In the case of the pointlike part ( |20| ) this bare distribution does, however, depend on 
the scale M, and its derivative with respect to InM 2 generates the inhomogeneous term in @. 
Formally g^g(x, Mo, M) as given in (18.2C) can be considered also for M < Mq, but is negative 
there. For M/Mq — > oo the second term in brackets of (18) vanishes and therefore all the pointlike 
solutions share the same large M behaviour 

q^(x,M ,M) - -^-a m (x) = q^(x,M)<x\n^-, (21) 

defining the so called asymptotic pointlike solution q^(x,M) ||, 0. Note that this asymptotic 
pointlike solution is a very special case of the general pointlike one fll8D, because for this solution 
the lower integration limit Mq in (pC]) has been identified with A, i.e. Mq = A! The fact that for the 
asymptotic pointlike solution ( |2l| ) a s (M) appears in the denominator of ( |2l| ) has been the source of 
claims (see, for instance, ||]) that q(x, M) = 0{l/a s ). This claim is wrong for two reasons. First, it 
is obviously invalid for those of the currently used parameterizations that distinguish the pointlike 
and hadronic components of photonic PDF. For instance, the widely used Schuler-Sjostrand sets 
SaSl and SaS2, which take Mq = 0.6 GeV and Mq = 2 GeV, respectively are not asymptotic 



pointlike solutions and therefore manifestly do not behave as a/a s . However, as argued in [14| this 
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Figure 2: Comparison, left in the nonsinglet and right in the singlet channels, of the functions 
xk(x) given by the splitting term (a/27r)3e 2 .(x 2 + (1 — x) 2 ), with the functions xa(x) corresponding 
to asymptotic pointlike solutions (|2TD. 



claim is misleading even for the asymptotic pointlike solution. To see the point consider the limit 
A — > for fixed M> We easily see that 

q^(n, M, M ) - ^%s(n, M) In (22) 



i.e. (18) reduces to moments of the first term in the expansion (P0|), corresponding to pure QED 
splitting 7 — > gg. The fact that the asymptotic pointlike solution (|2l|), for which Mo = A, diverges 
when A — > is then a direct consequence of the fact that for this (and only this) pointlike solution 
the decrease of the coupling a s (M/A) as A — > is overridden by the simultaneous extension of the 
integration region as Mo — > 0! In other words, if QCD is switched off by sending A —* without 
simultaneously extending the integration region, i.e. for fixed Mo, there is no trace of QCD left and 
we get back the simple QED formula. 

In summary, the pointlike part g^g(x, M) of the quark distribution function of the photon is of 
the order a and not of the order a/a s , as often claimed. One can use the latter only as a shorthand 
for the specification of large M behaviour as expressed in (^H). In fact, this is what one finds in 
the original papers |l^, |l3| ^ which do not contain any explicit claim that q(x, M) = 0(l/a s ). To 
take this behaviour literally leads, as shown in the next Section, to incorrect conclusions. 



3 Factorization scheme analysis of (x, Q 2 ) 

First a general comment. All the splitting and coefficient functions discussed below can be calcu- 
lated using one of two different approaches 

• Ultraviolet renormalization of composite operators within OPE. 
10 See, for instance, eq. (3.20) of (5J) 
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• Infrared renormalization technique anchored in the framework of parton model and based on 
analysis of Feynman diagrams [[!(]] . 

For the discussion of the transition between the properties of real and virtual photon, discussed 
in the next Section, the second approach has clear advantage in its simplicity and physical trans- 
parency and we shall therefore adopt it also throughout the paper. Moreover we shall work with 
nonzero quark masses m q , which provide natural regulators of parallel singularity associated with 
the primary 7 — ► qq splitting. A simple and straightforward analysis of this splitting (see next 
Section) yields 



C^{x,Q/M) 



(x z + (l-xY) I In -^ T + In — — I +8x(l-x)-l 



(23) 



Note that when evaluated via ultraviolet renormalization of composite operators using dimensional 
regularization contains also a term proportional to ln47r — 7^. However, this term has nothing 
to do with QCD and is related exclusively to the renormalization of electromagnetic couplant aP[ 
The reason why the existing NLO analyses of F^(x, Q 2 ), like |p], 15, 16j, |l7], |l8, IS], are incomplete 
is that they all take literally the relation q cx l/a s . This, in turn, leads to several incorrect 
conclusions: 



The term in (|12[) , which starts at the order 0(a), is claimed to be of the NLO in a s with 
respect to q(x,M), while in fact it is of the same order as q(x,M). 

The term in ( |i~5"l) proportional to aC^ is retained (though misleadingly assigned to NLO), 
while the term proportional to aa s C^ is discarded, while in complete NLO analysis it must 
be retained as well. 

Similarly in the case of splitting functions (PHi"l|) : while in the homogeneous splitting functions 



terms proportional to a 2 P^ are retained, the inhomogeneous splitting functions of the same 



order, a 2 kq 2 \ a 2 k^ are discarded, while they should be kept as well. 
As a result, all existing analyses write F^ at the LO in the form 

F^ LO (x,Q 2 ) = ^2xe^(x,M), (24) 
1 

where q[x, M) satisfy LO evolution equations, which take into account only kq° and pfj splitting 
functions and at the NLO use the expression 

F 2 ^ NLO (x, Q 2 ) = J2 2xe 2 q q{x, M)+ (25) 
q 

X; 2xe 2 [q(M) ® CW(Q/M) + G(M) ® c£\q/M)] + £e^f (x, Q/M) 

where the evolution equations for quark and gluon distribution functions include terms up to 
k?q\ k^Q and P, fj '. However, taking into account that q(x,M) = 0(a), correct LO and NLO 
expressions for read instead as follows 

F2' L °(x, Q 2 ) = X 2xe 2 q Ux, M) + ^e^f (x, Q/M)) , (26) 

11 Sec also discussion following eq. (^8|) below. 
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F^ NLO (x,Q 2 ) =£>ej Uaj,M) + ^effi\x,Q/M)\ + (27) 



where PDF in (^7j) satisfy evolution equations with splitting functions (pj-|rT|) up to the order 
a 2 (M), i.e. including terms proportional to kq,k(p. The difference of (|26|-[27|) with respect to 
dl-H) is threefold: 

• The term proportional to Oy is part of the LO expression. 

• The NLO expression in contains, beside the usual partonic convolutions q ® Cq 1 ^ and 
G (8> Cq also the NLO direct term proportional to C^~\ This additional term is crucial 
for the consistency of the NLO approximation. As shown below, it also cancels part of the 
factorization scale and scheme dependence of the NLO hadronic terms and thus contributes 
to theoretical stability of NLO calculations. 

(2) 

• In the NLO evolution equations inhomogeneous splitting functions k\ ,i = q, G of the order 
a 2 are retained. 

We note that in the case of the proton the expressions for F2(x,Q 2 ) at LO and NLO differ from 
those in ( ^6|]27| ) merely by the absence of the terms proportional to Oy and C^. For F^x^Q 2 ) 
factorization scale and scheme dependence of parton distribution functions is cancelled by explicit 
factorization scale and scheme dependence of the coefficient functions Cq , Cq ,i > 1. Concretely, 
at the NLO the scale and scheme dependence of quark distribution functions in the first (i.e. 
LO) term on the r.h.s. of (^) is cancelled to order a s by the corresponding dependence of NLO 
coefficient coefficients Cq 1 ^ , Cq^ . This reflects the fact that logarithmic derivatives of PDF of 
the proton with respect to InM and FS start at the order a s and also implies that at the LO the 
relation between F£(x, Q 2 ) and quark distribution functions is unique. Moreover, a precise physical 
interpretation of the factorization scale M is of little relevance because any redefinition of M can 
be compensated by appropriate change of the NLO splitting functions pff . 

For photon the inhomogeneous term in evolution equations for quark distribution functions 
modifies this cancellation mechanism because part of the factorization scale dependence driven by 
this inhomogeneous term is compensated already at the LO by the term in C^ proportional to 
kq ln(Q 2 /M 2 )! The presence of such a term in Cy , as well as the fact that c\^ contains the 
term proportional to kq , are general consequences of the RG invariance of F^ix, Q 2 )- For photon 
the interpretation of factorization scale M therefore does matter, as it determines LO coefficient 
C^\ Note that M 2 in the (|23| ) must be interpreted as maximal transverse momentum squared of 



the virtual quark included in the resummation (p0|). 

To see explicitly how the cancelation of factorization scheme dependence works at the NLO and 
why it requires the inclusion of the term C^ , let us follow the standard procedure Q of redefining 
quark distribution functions of the photon, but include terms up to the order aa s : 

q(x,M) = q(x,M) + ^e 2 Mx) + ^-^l e 2 q E,(x). (28) 



It is straightforward to show that provided q{x, M) satisfies the evolution equation with inhomo- 
geneous splitting functions kq ^,^ 1 ^ and kq 2 \ the tilded distribution function q(x,M) satisfies the 
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same evolution equation, but with tilded inhomogeneous splitting and coefficient functions given 
as: 



kf) 



feW- 



P 2p(0) 



7' 



e\ (p« ® D 7 + P« ® £ 7 + /? ^ 7 /2) . 



'7) 



(29) 
(30) 
(31) 
(32) 
(33) 



On the other hand, NLO homogeneous splitting functions P-p are unchanged by the redefinition 
(Pq) . The relations (|30|-31) imply that instead of the functions D 1 {x) and E 7 (x) this redefinition 

can equally well be parameterized by the inhomogeneous splitting functions k^p , . 

Despite the fact that D 1 is related by (|30|) to the splitting function k^ standing in @ by a s , 
only the term proportional to aa s E^ is related to genuine QCD effects. This is indicated already 
by the observation that solutions of the equation 

dq(x,M) 



a 



4°)( 3 



(34) 



dlnM 2 2vr 

to which the evolution equation (||) reduces in the limit a s — ► are determined up to an arbitrary 
function of x, i.e. just like the term in (|2q ) proportional to D^(x). The fact that this term in (pq) 
has nothing to do with QCD follows also from the fact that it can be fully included in the redefinition 



of the lower integration bound Mq in (|20|) by setting Mq(x) = Mq exp(—D 1 (x)/ (x 2 + (1 - x) )). 



With such a choice of the lower integration bound the result of the resummation in ( |20| ) does not 
vanish at any fixed initial scale, but still solves the same inhomogeneous evolution equation and 
therefore represents legitimate pointlike solution. 

It is also illustrative to see how the redefinition (p£|) affects separately the pointlike and hadronic 
parts of quark distribution functions. For the term containing the answer can be taken from 

the analysis in JlfJ, where explicit expressions for moments i^(n, M 2 ) as functions of kg(n) are 



given. One finds that the variation of k, 



(i) 



by 5k { q l) 



k 



(i) 



k, 



pointlike part q PL (n, M) at all scales M by the term (d(n) = 2Pqq' / @q 



(i) 
q 

(0) 



-p 2 P 



(0) 



modifies the 



5q PL (n,M) 



a 
2^ 



a s (M) 
a s (M ) 



-d(n) 



6k 



(i) 



a 
2^ 



a s (M) y d(n) 

(Mo); 



e g L> 7 (n), 



while in the hadronic part only the initial condition at Mq is changed 



5q HAD (n,M ) 



a 



5k 



(i) 



a 



2 *p£\n) 



—e 2 q D y (n). 



(35) 



(36) 



The NLO splitting functions remain, however, unchanged. Summing ( |35[ ) with fl3q) m ultiplied by 
the evolution factor (a s (M) / a s (Mo))~ d ^ yields the second term on the r.h.s. of Q28|). 



4 Structure of the virtual photon 

4.1 Equivalent photon approximation 

All present knowledge of the structure of the photon comes from experiments at the ep and e + e~ 
colliders, where incoming leptons act as sources of transverse and longitudinal virtual photons ^|. 
In a typical "photoproduction" experiment at HERA the average photon virtuality <P 2 ) = 10" 2 - 10" 3 GeV 2 . 
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q(pK Q ) 



H(Ph) 




b) 



c) 



Pt 



Figure 3: Kinematics of the basic j*(P 2 



) ~~ > QQ splitting in QED. 



To order a their respective unintegrated fluxes are given as 



a ( \ + {l-y) 2 ) 1 



2m 2 



y 



) 



(37) 



fl(y,P 2 ) 



2ir\ y P 2 

a 2(1 -y) 1 
2tt y P 2 ' 



p4 



(38) 



The transverse and longitudinal fluxes thus coincide for y = 0, while at y = 1, /2 vanishes. Note 
that while for P 2 ^> m 2 the second term in ( |37l) is negligible with respect to 1/P 2 , for P 2 close to 
^min = ™- 2 2/ 2 /(l — 2/) their ratio is finite and approaches 2(1 — y)/(l + (1 — y) 2 )- 

4.2 From virtual to real photons 

The choice of a way mass singularities resulting from the primary 7* — > qq splitting are regularized 
is crucial for smooth and physically transparent transition between the properties of virtual and real 
photons. As far as perturbative calculations of coefficient functions (hard scattering cross-sections 
in general) are concerned, there is no principal difference in this respect between QED and QCD. 
For PDF of the photon the situation is, however, more complicated. The point is that in QCD the 
transition from virtual to real photon is expected to be determined by nonperturbative parameters 
related to colour confinement, rather than the quark masses, which govern such transition in pure 
QED. Nevertheless, it is illustrative to see how this transition is realized in QED, where the "lepton 
distribution functions" of the photon, are explicitly calculable. 

In QED masses of charged fermions, in particular electron, play a fundamental role. For real 
photon parallel singularity associated with the purely QED splitting 7 — > qq is naturally regulated 
by quark and lepton masses. For virtual photon also its virtuality P 2 shields off this singularity 
but quark and lepton masses are indispensable for proper limiting behaviour of lepton (quark) 
distribution functions of the virtual photon as P 2 — ► 0. The effects of nonzero photon virtuality 
P 2 are threefold: 

• longitudinally polarized photons must be taken into account, 

• unintegrated PDF of transversally polarized photons contain terms proportional to P 2 , 

• parton level cross-sections obtain contributions from terms proportional to P 2 . 
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4.3 Parton model interpretation of constant terms in Cj^' 

The evaluation of quark distribution functions of the virtual photon in pure QED serves as a useful 
guide to parton model interpretation of the so called "constant terms" in LO coefficient function 
C-y and leads to explicit formulae for the virtuality dependence of quark distribution functions of 
the photon. 

Consider, for definiteness, the production of a heavy scalar particle H with mass Mjj in electron- 
proton collisions. The relevant parton level hard scattering subprocess 

j(P 2 ,k a )+q(p) -^e(k c )+H(p H ) (39) 

is described by the sum of amplitudes corresponding to diagrams in Fig. 3a, b. In pure QED and to 
order a the probability of finding inside the photon of virtuality P 2 a quark with mass m q , electric 
charge e q , momentum fraction x and virtuality r = m 2 — k 2 up to M 2 is given by the t-channel 
diagram in Fig. 3a 



n \ r M WAx m 2 P 2 ) 
g QED (*,m 2 ,P 2 ,M 2 ) = ( |-3e 2 ) J ^ Wt( ffi > dr ( 10, 



In the collinear kinematics, which is relevant for finding the lower limit on r, the values of m q , x, r 
and pt are related to initial photon virtuality P 2 as follows 

2 2 

1 — x 1 — X 

In terms of Mh and s = (k a + p) 2 , x = Mfj/s. W(x, m 2 , P 2 ,M 2 ) can in general be written as 

2 2 

W(x,m 2 P 2 ,M 2 ) = f(x)^—+g{x)m 2 a + h(x)P 2 + c(x)— + --- 
H 1 — X H s 



f(x)T + (g(x) - ^) m\ + (h(x) + xf(x)) P 2 + c(x)^ + • • • (42) 



where the dots indicate term of the type r k+1 /s k ,k > 1. There is no term proportional to s as 
it would violate unitarity. The functions f(x),g(x) and h(x) are unique functions that can be 
determined from the analyses of the vertex 7 — > qq in collinear kinematics. On the other hand, c(x) 
is a process dependent. The terms in ([42] ) proportional to f(x),g(x) and h(x) are dominated by 
small quark virtualities and have therefore clear parton model interpretation: so long as r <C Mfj 
eq. ([!(]) describes the flux of quarks that are almost collinear with the incoming photon and "live" 
long with respect to the production time of the heavy particle H in the lower vertex in Fig. 3a. 
Substituting ( f42| ) into ( f40| ) and performing the integration gives, in units of 3e 2 a/27r, 



M 2 \ 

q Q ED(x,m 2 q ,P 2 ,M 2 ) = /(x)ln — - + 



g(x)m 2 + h(x)P 2 
-f{%) + — 



1 




^-111111 

(43) 

In practical applications the factorization scale M is identified with some kinematical variable 
characterising hardness of the interaction, like Q 2 in DIS or Mjj in our process (39). As a result 
M 2 1 s becomes a function of x and we can write c(x)M 2 /s = k(x). For r mm = xP 2 + m 2 / (1 — x) -C 
M 2 (H) simplifies to 



9 „9 „,?n p / \ i f M 2 \ , g(x)m 2 + h(x)P 2 

wed(». < ' ) - /(*) '° xP> + mg/(1 _ I} - / w + + mV ^ 1 x) + - w- (i4) 



The s-channel diagram in Fig. 3b is, however, crucial for preserving gauge invariance in the case of the longitudinal 
photon where it cancels one of the terms originating from Fig. 3a that does not vanish for P 2 — > 0. 
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This basic result, which takes into account nonzero quark mass m q as well as initial photon virtuality 
P 2 , reduces for x(l — x)P 2 3> m 2 to 



q QED (x,0,P 2 ,M 2 ) =/(x)ln 



M 2 
xP 2 



el \ H X ) I \ 

f(x) + -^- + K(x), 



(45) 



whereas for P 2 /m 2 — > it approaches 



qQEl)(x,m 2 0,M 2 ) = /(x)ln 



'M 2 (l 



m„ 



f(x) +g(x)(l - x) +k(x), 



(46) 



describing distribution of quarks inside the real photon. The fact that (^) requires nonzero quark 
mass reflects the fact that in QED masses of charged fermions play vital role and cannot be zero. 

As in the case of the photon fluxes (37-3^) the leading logarithmic term, dominant for large 
M 2 , as well as the "constant" terms proportional to f(x),g(x) and h(x) come entirely from the 
integration region close to r mm and are therefore unique. At r = r mm both types of the singular 



terms in (|40|), i.e. 1/r and m 2 /T or P /t , are of the same order but the faster fall-off of the 1/t 
terms implies that for large M 2 the integral in ( |40| ) is dominated by the weaker singularity 1/r. In 
other words, while the logarithmic term is dominant at large M 2 , the constant terms resulting from 
nonzero m 2 and/or P 2 come from the kinematical configurations which are even more collinear than 
those giving the logarithmic term. On the other hand, not all constant terms are of this origin, 
as exemplified by the term proportional to k(x) which comes from the integration over the whole 
phase space and is therefore process dependent. These terms have no parton model interpretation 
and belong therefore naturally to the coefficient function Cy . 

The analysis of the vertex ^(P 2 ) — ► qq in collinear kinematics or the explicit evaluation of the 
diagrams in Fig. yields the following results 



Mx) 



x 2 + 



0, 



9t(x) 
9l{x) 



i 

l-x ' 

0, 



h T (x) 
h L (x) 



0, 

4x 2 (l 



x). 



(47) 



The vanishing of Jl(x) and <?l(x) is a consequence of gauge invariance, which guarantees that jl 
decouples in the limit P 2 — > 0. The fact that hx{x) = is due to helicity conservation, which 
permits violation by quark mass terms only. The expressions (42-44) exhibit explicitly the smooth 
transition between the quark distribution functions of the virtual and real photon in pure QED, 
governed by the ratio P 2 /m 2 . 

For virtual photon eq. (|l2;) holds separately for both transverse and longitudinal polarizations. 
For x(l — x)P 2 S> m 2 the finite terms CT, Clf are given as 



C7^(x,P 2 ,g/M) 
C[%(X,P 2 ,Q/M) 



(x 2 + (1 - x) 2 ) In + In + 8x(l 



Q- 



x) 



Ax{l - x) 



while for P 2 = 



C®(x,0,Q/M) 
d°|(x,0,Q/M) 



(x 2 + (1 _ x)2) L^ +]n l_£ 



+ 8x(l -x) - 1 



(48) 
(49) 

(50) 
(51) 



In combination with ([i4|) and ( |47| ) these results imply k(x) = —1 + 6x(l — x). The origins of the 
nonlogarithmic parts of in (48-50|) can then be identified as follows: 
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Real photon: 

-l + 8x(l-x) = -2 + 8x(l - x) + 1 
«■ » ' 

for massless quark 9 ( x )(^-^ x ) 

= -1 + 6x(l - x) - (x 2 + (1 - x) 2 ) + 1 (52) 

* v ' " W ' 

nonuniversal part f T ( x ) ^ { x )0-~ x ) 

Virtual photon: 

- 2 + 8x(l - x) = -1 + 6s(l - x) - (x 2 + (1 - x) 2 ) (53) 

V v ' " v ' 

nonuniversal part f T ( x ) 

As the virtuality P 2 (or more precisely r mm ) increases toward the factorization scale M 2 , the 
expression (|43|) for the quark distribution functions of the virtual photon approaches zero. We 
emphasize that this holds not only for the logarithmic term but for all terms that have parton 
model interpretation. 

4.4 The real world 

In realistic QCD nonperturbative effects, in particular those connected with the confinement, are 
expected to determine the structure of the real photon as well as its virtuality dependence. For 



instance, within the Schuler-Sjostrand set of parameterizations [21] the role of such parameter is 
played by vector meson masses for the hadronic component and by the initial Mq for the pointlike 
one. Nevertheless, the analysis of the previous subsections is still relevant for the coefficient function 
Cy as well as the discussion of virtuality dependence of the pointlike part of quark distribution 
functions. 

5 Phenomenological implications 

The results of previous Sections have several implications for phenomenological analyses of hard 
scattering processes with photons in the initial state. Some of them are discussed below. 

5.1 LO and NLO analysis of F^x^Q 2 ) for the real photon 

As argued in detail in Section 3, the existing analyses of F%{x, Q 2 ) are incomplete due to the fact 
that they do not take into account 



the NLO photonic coefficient function C^\x) in expression (27) for an 



• the NLO inhomogeneous splitting functions kq^ and Kq in evolution equations (|3[]5|). 

There is no problem to remedy the first shortcoming, because the coefficient function C-y (x) is 
actually known. As argued in [^] it can be deduced from the 0(a 2 ) gluonic coefficient function cal- 
culated in [23]. For illustration of its numerical importance, (a s (M)/vr)cf } (x) in MS+MS scheme 



is compared in Fig. |4|a to 2C^°^(x) as well as to (4-7T / ae 2 )u{x) and (4ir / ae 2 )u(x) corresponding to 
SaSIM parameterization. In this comparison we took a s (M)/2TT = 0.03, which is approximately 
the value of a s at M = 10 GeV. Note that the contribution of to F^{x, Q 2 ) is comparable to 
that of Cy . As in most experiments Q 2 is well below 100 GeV 2 , the relative importance of C^ 1 is 
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Figure 4: Comparison of the contributions of coefficient functions Cj (x) and C^' (x) to i 7 ^ with 
those of quark distribution functions u{x) and d(x) (a). In b) the SaSIM parameterization of xu(x) 
(solid line) is compared to the quantity fl54[) (dotted curve). 



even bigger. We also see that for x ~ 0.4 sum of the contributions of coefficient functions Cj and 
Cy to F2 is almost the same as the contributions of quark distribution functions themselves. 



As far as the NLO inhomogeneous splitting functions kq (x) and k^' (x) are concerned, the 
situation is different. These functions are unknown as they cannot be derived from the existing 
NLO calculations. To obtain them requires complete NNLO calculation, similar to that in [24] for 
moments of nucleon structure functions. No such calculation is, however, in sight. In the absence 
^'x) and Kq (x) a complete NLO analysis of i 7 ^ is impossible, but we may investigate what 



.(2), 



of k 



happens if at least the effects of C 7 
effects is to evaluate the difference 



(i) 



are taken into account. The simplest way to estimate these 



q'(x,M) = q(x,M) 



a a s (M) 



2vr 2vr 

defining quark distribution functions which when inserted into 



(54) 

|]) including the NLO coefficient 

function C^' return the same K7 as q(x, M) inserted into (^7j) without the term. For u- 
quark the effect of the redefinition (p3|) is shown in Fig. Qb, using SaSIM parameterization of 
u(x,M) at M = 10 GeV. The reduction of u(x,M) due to the inclusion of Cj is substantial for 
all 0.05 < x < 0.9. 

In the situation when full NLO calculation of F2 is impossible due to the lack of knowledge 



of k q and kg' and Fig. ^ indicates the importance of C^' term, LO analyses of and other 
quantities involving initial photons is all we can consistently do. In such analyses we can exploit 
the freedom in the choice of the factorization scale and scheme. Contrary to the case of hadrons, 
these choices influence the relation between PDF and F^ (and other physical quantities) already 
at the LO through factorization scale and scheme dependence of . It is therefore not true, as 



.(2) 



(1) 



claimed for instance in [25], that the SaSIM and SaS2M sets of parameterizations are "theoretically 



inconsistent" because they combine in LO expression for F 2 7 the "NLO" quantity with the LO 
quark distribution functions. The MS sets of SaS parameterizations are as legitimate definitions 
of the LO quark distribution functions as their DIS companions and there is no theoretical reason 
why, as suggested in [29], they should not be used in phenomenological analyses. 
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5.2 What is measured in DIS on virtual photons? 

In experiments [||, ^ || || at e + e~ colliders structure of the photon was investigated via standard 
DIS on the photon with small but nonzero virtuality P 2 . The resulting data were used in ]l8| . pi 
to determine PDF of the virtual photon in LO and NLO approximations. In these analyses C 7 
was taken in the form 



C^\x,P\Q/M) 



M 2 1 \ 
+ (l-*) 2 )lln^+ln-J+6z(l 



(55) 



which, however, does not correspond to the structure function that is really measured in e + e 
collisions, but to the following combination 



1 



(56) 



of structure functions corresponding to transverse and longitudinal polarizations of the target pho- 
ton. This combination results after averaging over the target photon polarizations by means of 
contraction with the tensor —g^ v /2. The term —2 + 6x(l — x) follows also directly from consider- 
ations of the previous Section: 



2 + 6x(l-s) 



-2 + 8x(l -x)-2x(l -x) 

from 7 T from "f L /2 



(57) 



Because the fluxes (37]-38|) of transverse and longitudinal photons are different functions of y, 
any complete analysis of experimental data in terms of the structure functions K7 T (x, P 2 , Q 2 ) 
and F^ L (x,P 2 ,Q 2 ) at fixed x,P 2 ,Q 2 requires combining data for different y. This is in prin- 
ciple possible, but experimentally difficult to accomplish. The situation is simpler at small y, 
where fj^(y,P 2 ) = f2(y,P 2 ) = f^(y,P 2 ), and the data therefore correspond to the convolution 
of f^(y,P 2 ) with the sum F^ T + F^ L - ^he nonlogarithmic term in appropriate to this com- 
bination is, however, not —2 + 6x(l — x), as in (55) and [|l^, 19] but —2 + 12x(l — x), the sum 
of nonlogarithmic terms corresponding to transverse and longitudinal photons. Numerically the 
difference between these two expressions is quite important. 



5.3 Longitudinal gluons inside hadrons? 

If there are longitudinal photons inside leptons, what about longitudinal gluons inside hadrons? 
Although real longitudinal gluons decouple, as do real longitudinal photons, gluons as well as quarks 
inside hadrons are off-shell and therefore there in no reason not to introduce distribution functions 
of transverse and longitudinal gluons as well. The problem is, however, that while in QED the fluxes 
fj<(x, P 2 ) and fi(x, P 2 ) are known, we do not know how to calculate analogous fluxes of transverse 
and longitudinal gluons inside hadrons. In fact it would be sufficient to know the relative size of 
these fluxes, but even this is not calculable. Nevertheless, guided by the situation for photon fluxes, 
we can expect them to be similar, in particular at low x. If this is assumed, the considerations of 
the preceding subsection applied to gluons imply that the NLO gluonic contribution to F^x^Q 2 ) 
comes from the sum of contributions of transverse and longitudinal gluons. This in turn means 
that the finite nonlogarithmic term in the gluonic coefficient function Cq(x) should be taken as 
—2 + 12x(l — x) rather than the usual —2 + 6x(l — x). 



15 



5.4 Virtuality dependence of initial conditions 

There are currently two approaches to introducing virtuality dependence of PDF of the photon. 
Schuler-Sjdstrand base their parameterizations pl[ on the idea, suggested in [{26(1 , to use dispersion 
relations in P 2 written for moments of F^(x, P 2 , Q 2 )- Their parameterizations do not satisfy the 
same evolution equations (|3[|5|) as those of the real photon, but as the difference is formally of power 
correction type this is no principal drawback. 

In a different approach pursued in |l8|, [l5|, PDF of the virtual photon are assumed to satisfy 
the same evolution equations @]5|) as those of the real photon, and their virtuality dependence is 
introduced via the virtuality dependence of the initial conditions at Q 2 = P 2 [^]. This dependence 
is assumed to interpolate between the standard boundary conditions for the real photon, which 
in the GRV approach are defined at a very low scale Q 2 = fi 2 < 1 GeV 2 , and the predictions of 
perturbation theory 

D(x,P 2 ,P 2 ) =7 ] (P 2 )D N p(x,P 2 ) + [l-7 ] (P 2 )}D PT (x,P 2 ), D = q,q,G, (58) 
where P 2 = max(P 2 ,fi 2 ), r](P 2 ) = (1 + P 2 /m 2 )- 2 and 

dm^p 2 ) = o<p >2 <\ 2 > ( 59 ) 



~ 2 f 0, LO 

q PT (x, P ) = k (W/p) | 2e 2o_ ^ + (1 _ x?) ln ^ _ 2 + 6x(1 _ x) ] f NLO (60) 

G PT (x,P 2 ) = 0, (61) 



where k, f p , ji , f n (x, \i ) are given in |27], p8|. The above boundary conditions are assumed to be 



valid in the DIS 7 factorization scheme (see next subsection). The difference in (|60| ) between the 
form of boundary conditions on quark distribution functions at LO and NLO follows again from 
misassignment of the coefficient function Cy (x) to the NLO. As argued in previous sections this 
term is actually of the LO and therefore if taken into account in boundary conditions, it must be 
present in both LO and NLO analyses. 

To include it would, however, go against the very essence of the parton model interpretation 
of virtuality dependent PDF. To retain clear physical interpretation of the factorization scale Q, 
q(x, P 2 ,Q 2 ) must vanish as the lower bound on the quark virtuality r mm (x) oc P 2 approaches Q 2 
from below. This occurs explicitly in the expression (f43|) for all terms, whether logarithmic or not, 
that have parton model interpretation. The term r mm /Q 2 , which guarantees the vanishing of the 
nonlogarithmic part of ( |43| ) as r mm — > Q 2 , is formally of the higher twist and therefore invisible in 
||] or any other considerations within the leading twist approximation. 

5.5 Why DIS 7 factorization scheme? 

For hadrons the LO relation between quark distribution functions is the same in all FS and is 
also identical to that of the parton model. DIS FS was introduced to retain this simple relation 
also at the NLO, where different FS correspond to different NLO splitting functions pjp(x) and 
hard scattering cross-sections (coefficient functions for F2). However, once the latter is chosen, the 
former is uniquely determined (or vice versa). For proton the DIS FS amounts to setting = 0. 

For photon the situation is different. The DIS 7 FS was introduced for real photon to get rid of 
the troubling ln(l — x) term in the LO expression (f26|) for (x, 0, 1). It implies setting Z? 7 = C 7 °^ 
14 To retain the notation of [jl s|, |l9| I use in this subsection Q 2 instead of M 2 to denote the factorization scale. 
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Figure 5: Exact as well as approximate bounds on the ratio r/m 2 for real photons and massive 
quarks (a) and on the ratio t/P 2 for virtual photons coupled to massless quarks (b). In both cases 
solid and dashed curves correspond to exact and approximate bounds respectively, plotted for four 
values of the ratio Q 2 /m 2 in a) and Q 2 /P 2 in b). 



and Ery = in (28) and thus modifies the inhomogeneous splitting functions and k^ 2 \ but 

does not change the NLO homogeneous splitting functions Pm ■ As emphasized in Section 3 this 
redefinition of quark distribution functions has in fact little to do with QCD. In parton model 
derivation of Section 4, the troubling term ln(l — x) is a direct consequence of the fact that for 
real photon parallel singularity associated with the splitting 7 — * qq is regulated by means of 
quark masses. This is natural in QED, but as emphasized in subsection 4.2 not in QCD, where 
nonperturbative properties of hadrons (with the exception of pseudoscalar mesons) and the photon 
are expected to result from the dynamics of confinement and we can therefore start from massless 
quarks. 

Note that the term ln(l — x) in Cj(x, 0, 1) causes problems primarily because it appears there 
decoupled from the value of the quark mass m q with which it originally entered the considerations 
in the expression ( |4l| ) for r imn and thus persists there even in the limit m q — > 0. It comes from the 
lower bound on the quark virtuality r in collinear kinematics 



<T<T^=^-. (62) 



coll _ m q s ^ s noil _ Q 



111111 1 - X 

The ratio Q 2 (l — x)/xm 2 of upper an lower integration limits in (62) equals unity (and thus the 
integral vanishes) at x m = 1/(1 + m 2 /Q 2 ), which for light quarks is very close to 1. For x — ► 1 
collinear kinematics is no longer appropriate and should be replaced with the exact one 

y 1 - y/l-z t yl + VTTz 

2 x ^| - 2 x ' (63) 

where y = Q 2 /m 2 and z = 4yx/(l — x). These exact bounds are shown, for several values of y, 
in Fig. |5|a by the solid curves, together with the dashed ones, corresponding to the approximate 
bounds as given in (|62|). Also shown by the dotted curves are the functions 1/x and 1/(1 — x). 

On the other hand when the term ln((l — x)/x) alone is added (as part of the expression for 
C^) to the pointlike solution of the inhomogeneous evolution equation with initial conditions at 
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Mo 3> rriq, the sum vanishes already at xm = 1/(1 + Mq/Q ) <C £ m an d becomes independent of 

m g down to m g = 0! However, it would be strange to retain in Cy the term that follows directly 
from assuming nonzero quark masses, but which when combined with the PDF is then unrelated 
to their values. In the presence of confinement it seems more appropriate to start from the very 
beginning with massless quarks coupled to off-shell photons and then construct the limit P 2 — ► 0. 
This approach is even more natural taking into account that all current data on the structure of 
the "real" photon actually come from interactions of virtual photons, albeit with small virtuality. 
For them the bounds on quark virtuality in collinear kinematics 

O 2 

xP 2 <t <— (64) 

x 

are only slightly changed using exact kinematics 



(l-x)y(l- VT+z) t (1-^(1 + ^/1+1) 

2 x ~ P 2 ~ 2 x 1 ' 



where now y = Q 2 /P 2 and z = 1/(1 — x) — xy/(l — x). The modification of ( pip is so small that 
some of the solid curves in Fig. ||b are indistinguishable from the corresponding dashed one. Our 
suggestion is therefore to use also for the real photon the expression 

C( 0) (j;,0,1) = [rE 2 + (l-x) 2 jln4r-2 + 8x(l-x), (66) 
' L J x z 

derived for the transverse photon with virtuality P 2 S> m 2 . This term is different from (p5|) 
normally used for the virtual photon, as the latter corresponds to F^ T — F^ L /2. 



6 Conclusions 

We have discussed the factorization mechanism for Fjix, Q 2 ) and pointed out the differences with 
respect to the case of proton structure function F2(x,Q 2 ) which are due to the presence of the 
inhomogeneous term in evolution equations for g 7 . On the one hand, this term allows us to 
calculate the asymptotic behaviour of -F^(x,Q 2 ) as Q 2 — > oo but on the other hand the same 
term also implies that a complete NLO analysis of F^XjQ 2 ) requires the inclusion of the 0(a 2 ) 

(2) (2) 

terms k\ , k G in the inhomogeneous splitting functions k q , kc as well as of the 0(a s ) term 
term in photonic coefficient function Ci. Unfortunately, none of them has been included in 
the existing phenomenological analyses. There is no problem to include the latter as the necessary 
calculations are available, and we have shown that its numerical importance may be quite large. On 

(2) (2) 

the other hand, kq and k G are not known as their evaluation requires three loop QCD calculation. 
Consequently, at the present time a complete NLO analysis of F^{x, Q 2 ) is impossible to perform. 

As far as the structure of the virtual photon is concerned, we have discussed the question of 
what is actually measured in DIS on the virtual photon and emphasized the role of the longitudinal 
photon in these considerations. 

We have also analyzed parton model interpretation of the constant terms in C-^ and discussed 
their implications for the specification of initial conditions imposed on the PDF of the virtual 
photon. We have argued that the presence of colour confinement washes out the difference between 
the regularization of mass singularities in the case of real and virtual photons. This offers us a 



18 



simple way of avoiding the problems with the term ln(l — x) appearing in the standard expression 
for in the case of the real photon. 
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